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We consider a composite particle formed by two fermions or two bosons. We discover that 
composite behavior is deeply related to the quantum entanglement between the constituent particles. 
By analyzing the properties of creation and annihilation operators, we show that bosonic character 
emerges if the constituent particles become strongly entangled. Such a connection is demonstrated 
explicitly in a class of two-particle wave functions. 
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■ In Nature many particles are composite particles formed by two or more fundamental fermions. Atoms, molecules 
and nuclei are familiar examples of composite particles. As long as physical processes do not reveal the internal 

^ . structure, a composite particle as a whole can be treated approximately as a point-like boson (fermion) if the number 
O of constituent fermions is even (odd). Intuitively, the more compact the composite particle is, the better it behaves 
as a pure boson or fermion. This raises interesting questions how the degree of compositeness is quantified, and 
what physical effects would appear if the internal structure becomes apparent. In the case of particles formed by 
two fermions, several authors have addressed the compositeness effects of atoms in Bose-Einstein condensates 0, Q 

^ ^ and electron-hole pairs in semiconductor In particular, Combescot and Tanguy have indicated the deviation of 

^ ' purely bosonic behavior, based on the properties of the creation operator associated with composite particles 

, The issue of compositeness also appears in particles made of bosons. An important example of this kind is biphoton 
' states generated by optical spontaneous parametric down conversion. Such highly correlated photons have been 
, employed in fundamental experiments, e.g., to test the violation of Bell's inequality 0, and to realize Einstein- 
Podolsky- Rosen paradox [5|. Recent experiments have also explored the 'binding' effects revealed from the de Broglie 
wavelength of biphotons [a, ■ Interestingly, although the constituent photons are non-interacting and they can be 

■ spatially separated, the composite nature of biphotons can still be detected by suitable interference schemes. This 
suggests that the concept of compositeness is not limited to mechanically bounded particles in position space. 

In this paper we examine the composite representation of two-particle systems from the view point of quantum 
information. We show that quantum entanglement provides an understanding of the origin of composite behavior. 
As we shall see below, the degree of entanglement between constituent particles determines how close a composite 
particle behave as a pure boson. This implies an interesting picture that constituent particles are somehow bounded 
^ by quantum entanglement. Mechanical binding forces are not essential, they serve only as physical means to enforce 
quantum correlations. Since the constituent particles can be correlated in many different ways, the composite repre- 
^ ' sentation is not limited to position space or momentum space. In this paper we will indicate the fundamental role 
of entanglement, based on the properties of creation and annihilation operators associated with composite particles. 
In contrast to previous studies that concern only fermionic constituent particles, we provide a general analysis that 
applies to bosonic constituents as well. 

To begin with, we consider a composite particle C formed by two distinguishable particles A and B. Both A and B 
are either fermions or bosons. Let '^{xatXb) be the wave function of the two-particle system, we can always express 
the wave function in the Schmidt decomposition form: 

oc 

^{xA,XB) = Y.VKc^i^\xA)4>'-n\xB) (1) 
n=0 

where the Schmidt modes 4>[^'' {(f'lf ^) form a complete and orthonormal set for particle A [B). Specifically, 

and A„ are defined by the eigenvectors and eigenvalues of the reduced density matrix of particle A, and can be 
obtained similarly from the reduced density matrix of particle B. Eq. (1) reveals the quantum correlation by showing 

the pairing structure explicitly. If the particle A appears in the mode (pi^^ , then with certainty the particle B must 
be in the mode (j)n ■ The distribution of A„ provides a measure of entanglement. This is usually discussed in terms 
of the entanglement entropy E = — A„ logj A„. However, a more transparent measure of entanglement is to count 
the 'average' number of Schmidt modes actively involved. The Schmidt number K. provides this information 

oo 

JC^l/J2^l (2) 

n=0 
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The larger the value of K, is, the higher the entanglement. In fact, K, equals the inverse of the purity of single-particle 
density matrix, and it also equals the linear entropy apart from a constant. A disentangled (product) state corresponds 
to /C = 1, i.e., only one term in the Schmidt decomposition. If a Schmidt decomposition contains M terms of equal 
weight (i.e., A„ = \/M for n = 1, 2, .., M), then we have IC = M which is exactly the number of mode pairs involved. 

In the second quantized representation, Eq. (1) corresponds to a state generated by a creation operator acting 
on the vacuum. Such a creation operator is defined by, 



= E VA^at 6t (3) 



ri=0 

where aj^ {blJ is the creation operator of the particle A [B) in the Schmidt mode (p^f'^ {(j)^^). The commutation 
relation between c and is given by, 

[c,ct]^l + sA, (4) 

where 5 — +1 if A and B are bosons, s = —1 \i A and B arc fcrmions. The operator A in Eq. (4) is defined by, 

00 

A = E ^" (4an + bibn) (5) 

n=0 

which has non-zero matrix elements, depending on the states involved. Therefore c and are not strictly bosonic 
operators. 

To examine the properties of c and c\ we consider a system containing two or more C particles. Let us define the 
iV— particle states by, 

'^^^^^ 7m'°^ (6) 

where xn is a normalization constant such that {N\N) = 1. In order to test how good the operator c behaves as 
a bosonic annihilation operator, we need to determine how the operator acts on |7V). This is given by a general 
equation: 



c\N)=aN^N\N -l) + \eN), (7) 

where is a constant, and the correction term \eiq) is orthogonal to \N — 1). Such a correction term is necessary 
because the set of |A^) states (with all N) is only a subset of the entire Hilbert space associated with the constituent 
particles. A non- ideal bosonic operator would inevitably cause transitions into the states not described by (6). 
From Eq. (7), the annihilation operator c is bosonic if the following two conditions are satisfied: 

aN ^ 1 (8) 
{eN\eN) ^ 0. (9) 

After some calculations for iV > 1, we have 



aiv = W^^ (10) 

V XN-1 

{eN\eN) = l-N^^ + {N-l)^^^. (11) 

Xn-1 Xn 

Therefore the conditions (8) and (9) are controlled by the ratio of normalization constants. An ideal composite boson 
emerges in the limit Xn±iIxn 1- 
For convenience, let us write 

_ J A,B are fermions , „^ 

~ \ are bosons. ^ ^ 

By carefully counting the states allowed for ferimons and bosons, it can be shown that Xn Xn given by |llj |. 

X^ = N\ E Ap,Ap,....Ap„ (13) 

Pn>Pjv-i>---->P2>Pi 

Xn = NI J2 Ap,Ap,....Ap„. (14) 

Pn>Pn-i>---->P2>Pi 
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The calculation of these summations can be quite complicated in general. In the case particles A and B are fermions, 
can be analyzed by methods discussed in Ref. p|. However, exact analytical expressions of x^v for both fermions 
and bosons in closed forms are rare. 

In this paper we consider a realistic class of two-particle wave functions that allows exact closed form expressions 
of xn- Such a class of two-particle wave functions is specified by the Schmidt eigenvalues, 

A„ = (l-z)z" n = 0,l,2... (15) 

where the parameter z is defined in the range < z < 1. This parameter determines how rapid A„ decreases with n. 
A representative example of this class of wave functions is (double) gaussian given by: 

«'(a;A,2:B) =AAe-(^^+^«)'/'^'e(^-^~^«)'/'^'. (16) 

Here A/" is a normalization constant, and ctc and are widths along the xa -\-xb and xa — xb directions respectively. 
The Schmidt decomposition of Eq. (16) gives A„'s in the form of Eq. (15) with 

^=f^^V- (17) 



The corresponding Schmidt modes in this example are simply the eigenfunctions of a harmonic oscillator. 

To evaluate x^ and x^: we let: pi = q^, P2 ^ Qn + Qn-i, Pa = Qn + Qn-i + qN-2,-- Pn ^ Qn + Qn-i + ■■■ + qi, 
then Eq. (13) and (14) become, 

oo oc oo oo 

X^ - iV!(l -Z)^Y.Y1 E E ^9i+2?2+3g3 + ...+JV,„ (^8^ 

91=0 52=0 ,jjv-l=0 gjv=0 
OO OO oo OO 

xf, = m{i -z)^Y.Yl E (19) 

91=1 «2=1 g„_i=lgjy=0 

These summations can now be carried out easily, 

(l-z)(l-z2)....(i_^JV) ^^^) 



iV!z^(^-i)/2(l-z)^ 

(] 

Hence the normalization ratios are given by. 



(1-Z)(l-Z2)....(l-Z^)- ^''^ 



xf+1 ^ (A^+l)(l-z) 

(1-z^+i) 
Xn+1 _ n{N + 1){1-z) 



(22) 



F ~ ^ /I N+^ \ ■ (2^^) 
Xn (1-z^+i) 

From Eq. (7) and (10), the normalization ratio determines the modification of Bose enhancement factor. The results 
Eq. (22) and (23) indicate that X;^-i-i/Xn > 1 ^^^^ Xw+i/x^w < 1- The difference between the two types of constituents 
can be understood, because bosons tend to stay together in the same state, while fermions do the opposite due to 
the exclusion principle. However, quantum statistics of constituent particles becomes less important as z approaches 
one, where the composite particle become a pure boson. In the example of double gaussian wave function (16), the 
z 1 limit correspond to the cases Gr ^ CTc or (Tc 3> cr^. We remark that that the particle number iV is a key factor 
in Eq. (22) and (23). A larger number of particle number would require the value of z to be closer to one in order to 
maintain the normalization ratio. 

Now we can make an explicit connection with quantum entanglement. For the Schmidt eigenvalues given by (15), 
the Schmidt number K, defined in Eq. (2) takes the form: 

l-z2 

(1 - zY 

which is a monotonic increasing function in the range < z < 1. By expressing z in terms of /C, both Xw+i/Xw 
and Xat+i/Xw '^^^ t)^ directly related to the degree of quantum entanglement. As JC increases, we find that both 
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X^+i/x^ s-nd X^+i/x^ approach one monotonically for the entire range of /C. In particular, it can be shown that 
for /C> TV, 

XN+i/XN « 1 + sN/JC (25) 

where s defined in Eq. (4). Since the value of JC corresponds to an effective number of Schmidt modes |^, bosonic 
particle description is valid when the effective number of Schmidt modes is much greater than the total number of 
composite particles. 

The discussion above indicates that composite behavior strictly depends on the degree of quantum entanglement, 
at least for the class of two-particle wave functions with Schmidt eigenvalues specified by Eq. (15). We note that 
Eq. (15) corresponds to a wide variety of wave functions, because Schmidt modes can be chosen from any discrete, 
complete and orthogonal set of basis functions. Indeed, Schmidt modes do not play any role in determining xn ■ It is 
the distribution of A„ that controls the compositeness of the particle. As long as the Schmidt eigenvalues are specified 
by (15), different forms of Schmidt modes correspond to the same degree of compositeness. This also shares the same 
feature with quantum entanglement. 

To summarize, our work here examine the foundation of composite representation of two-particle systems. Within 
the class of two-particle wave functions, we show that the origin of composite behavior is ultimately related to quantum 
entanglement between constituent particles. Therefore composite representation can be applied to strongly entangled 
particles, which are not limited to mechanically bounded systems. Finally, we remark that the wave functions 
satisfying Eq. (15) are not completely general, but owing to flexibility of the parameter z, these functions may be 
used to approximate wave functions in actual situations. This suggests that the relation between composite behavior 
and the degree of quantum entanglement can be quite general. However, a full analysis would require the study of 
Eq. (13) and (14) for arbitrary distribution of Schmidt eigenvalues, which is a topic open for future investigations. 



Acknowledgments 

The author thanks Prof. M.-C. Chu for discussions. This work is supported in part by the Research Grants Council 
of the Hong Kong Special Administrative Region, China (Project No. 400504 and Project No. 423701). 



[1] S. S. Avancini, J. R. Marinelli, and G. Krein, J. Phys. A. 36, 9045 (2003). 

[2] S. Rombouts, D. Van Neck, K. Peirs, and L. PoUet, Mod. Phys. Lett. A 17, 1899 (2002); S. Rombouts, D. Van Neck, K. 

Peirs, and L. Pollet, Europhys. Lett. 63, 785 (2003). 
[3] M. Combescot and C. Tanguy, Europhys. Lett. 55, 390 (2001); M. Combescot, X. Leyronas, and C. Tanguy, Eur. Phys. J. 

B. 31, 17 (2003); M. Combescot and C. Tanguy, Europhys. Lett. 63, 787 (2003). 

[4] Z. Y. On and L. Mandel, Phys. Rev. Lett. 61, 50 (1988); Y. H. Shih and C. O. Alley, Phys. Rev. Lett. 61, 2921 (1988). 
[5] J. C. Howell, R. S. Bennink, S. J. Bentley, and R. W. Boyd, Phys. Rev. Lett. 92, 210403 (2004). 
[6] J. Jacobson, G. Bjork, I. Chuang, and Y. Yamamoto, Phys. Rev. Lett. 74, 4835 (1995). 

[7] E. J. S. Fonseca, C. H. Monken, and S. Padua, Phys. Rev. Lett. 82, 2868 (1999); E. J. S. Fonseca, Z. Paulini, P. Nussenzveig, 

C. H. Monken, and S. Padua, Phys. Rev. A 63, 043819 (2001); K. Edamatsu, R. Shimizu, and T. Itoh, Phys. Rev. Lett. 
89, 213601 (2002). 

[8] R. Grobe, K. Rzaiewski and J.H. Eberly, J. Phys. B 27, L503 (1994). 

[9] C. K. Law and J. H. Eberly, Phys. Rev. Lett. 92, 127903 (2004). 
[10] The derivation of Eq. (11) requires the expectation value (Af|A|Ai'), which can be obtained by using the algebraic method 

in Ref. [3]. 
[11] In deriving Eq. (13), we write 

c^^ |0> = ^ -yA^^^T^J^TTTA^f (fci, fe, fcjv) |fci, fc2, fejv) 

where \kj) denotes the occupation with an A particle and a B particle in the Schmidt mode kj, and F{ki, k2, fcjv) is the 
weight factor for the state \ki, k2, fcjv). If these fe's have d same terms (e.g. k\ = k2 = ... = fc^), and all others terms are 
distinct, then the weigh factor should be N\/d\ x d! = N\. The d! in the nominator comes from the fact that contains 
operators a\b\. This argument applies to general sequence with any combinations of degenerate terms. 



